On the computational complexity of the Riemann mapping 



I. Binder* M. Braverman^ M. Yampolsky * 

Dept. of Mathematics Dept. of Computer Science Dept. of Mathematics 

University of Toronto University of Toronto University of Toronto 

February 1, 2008 



Abstract 

In this paper we consider the computational complexity of uniformizing a domain with 
a given computable boundary. We give nontrivial upper and lower bounds in two settings: 
when the approximation of the boundary is given either as a list of pixels, or by a Turing 
Machine. 



1 Introduction 

1.1 Foreword. 

Computational conformal mapping is prominently featured in problems of applied analysis 
and mathematical physics, as well as in engineering disciplines, such as image processing. In 
this paper we address the theoretical foundations of numerically approximating the confor- 
mal mapping between two planar domains. We obtain a lower bound on the computational 
complexity of an algorithm solving this problem, and show that this bound is almost sharp. 
To achieve the latter, we present a very space-efficient probabilistic algorithm for constructing 
such a mapping. 
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1.2 Background in computational complexity theory 

We present here some basic definitions and results from the computational complexity theory. 
A more comprehensive discussion can be found in [Sip, Papa]. 

The primary goal of the computational complexity theory is to classify different computa- 
tional problems into complexity classes according to their computational hardness. The basic 
abstract object here is a Turing Machine which for most purposes can be thought of as a 
program in any programming language. 

The complexity class P includes problems that are computable in time polynomial in the 
length of the input. Those are thought of as the "relatively easy" problems. Examples of 
problems in P include arithmetic operations, finding a shortest path in a graph and primality 
testing. "Difficult" problems, such as factoring integers or computing the optimal strategy for 
playing "Go" on an n x n board, are generally thought not to be in P. Whether a problem 
is in P or not is usually a good criterion in assessing its true hardness. By an analogy with 
P one can define the class EXP of problems solvable in time 2™ c for some c on input length 
n. Factoring integers is in EXP via an obvious exhaustive search. Playing "Go" optimally is 
also in EXP, since we can easily enumerate all possible games and compute an optimal path 
in time 2°( n \ Using a diagonalization argument, it is not hard to see that P C EXP (see 
e.g. [Sip]). 

The complexity class NP contains problems that are easy to verify, but may be hard to 
guess. More precisely, a predicate Q(x) is in NP, if there is a poly-time computable predicate 
R(x,y), where y has length polynomial in the length of x, such that Q{x) = By R(x,y). By 
an exhaustive search for y one sees that NP C EXP. There is a subclass of NP called the 
NP-complete problems, or NPC. Problems in NPC have the property of being the "hardest" 
in NP: if one could solve any problem in NPC in polynomial time, then one could solve all 
NP problems in polynomial time. 

One of the most famous NP-complete problems is the satisfiability problem SAT. The 
problem is the following: given a propositional formula (p(y) does it have a truth assignment 
yT such that 0(?/t) = 1- An example of a problem in NP that is thought to be hard but not 
NP-complete is the following. Given a pair of numbers m < n, determine whether n has a 
divisor between 2 and m. This problem can be used to factor integers. It is in NP since it can 
be formulated as 3k (1 < k < m) A k \ n. It is one of the Clay $1,000,000 questions of whether 
P = NP. 

A bigger class of problems is the class #P. It is the class of problems which are equivalent to 
counting the number of satisfying assignments for a given propositional formula - this natural 
complete problem for this class is denoted by #SAT. Obviously NP C #P, since to solve 
SAT we only need to know whether the number of its satisfying assignments is bigger than 
or not, which is easier than actually determining this number. 

The next class of problems is the class PSPACE - the class of problems solvable in space 
polynomial in the input size. It is easy to see that all the classes mentioned above are in 
PSPACE. On the other hand, PSPACE C EXP, since a machine with p(n) memory bits 
can have at most 2 p (") different configurations, and can run for at most 2 p(n ) steps without 
getting into an infinite loop. 



2 



The class of problems solvable in logarithmic space is a class of problems that are solvable 
in space O(logn) for input size n. Here the input and the output are read-only and write-only 
respectively. This class is denoted by L. By the same reasoning as PSPACE C EXP, we 
have L C P. A randomized version of L are the problems that can be solved correctly with 
error probability < 1/n in space O(logn) and time poly(n). This class is called BPL. 

Overall, we have the following chain of inclusions: 

L C BPL CPCNPC#PC PSPACE C EXP. 

By diagonalization, BPL 7^ PSPACE, and P 7^ EXP. No other separations are known. 

In recent years, some progress has been made in derandomizing BPL, that is in showing 
that there is a deterministic algorithm that does not require a lot of additional computational 
resources. We will need the following recent result on derandomization: 

Theorem 1.1. [Nis] There exists a deterministic algorithm for the following problem: 

Input: An n x n transition probability matrix M , an integer t, and a rational e. 

Output: A matrix A such that \ \A — M*|| < e. 

The algorithm runs in time poly(N) and space 0(log 2 N), where N = n + t + e _1 . 

We will also use a circuit complexity class. A circuit consist of inputs, logical gates, and 
an output. The gates are usually NOT (one input, one output), AND, and OR. The latter 
gates can either have two, or unboundedly many inputs. In the discussion below, any number 
of inputs is allowed. The size of a circuit is the number of gates used. The depth of a circuit 
is the number of gates on the longest path from an input to the output. It is known that any 
boolean function / : {0, l} n — ► {0, 1} can be computed by a circuit of size 0(2 n /n). It is not 
hard to see that functions in P are computable by polynomial size circuits. The class AC is 
the class of functions that are computable by a family of circuits (one for each input size) , that 
have constant depth and polynomial size. This is one of the very few complexity classes for 
which non-diagonalization lower bounds exist. In particular it has been shown that computing 
the parity of the number of l's in a string cannot be done in AC (see, for example, [FSS]). 
A more general problem that cannot be done in AC is the majority problem MAJ n : given a 
string x £ {0, l} n , MAJ ra (x) is 1 if and only if the majority of the entries in x are 1. 

1.3 Computational complexity of sets 

We review the definition and the basic properties of computable sets. We refer the reader to 
[BW, Wei, RW, Brav] for a more comprehensive exposition. 

Intuitively, we say the time complexity of a set S is t(n) if it takes time t(n) to decide 
whether to draw a pixel of size 2~ n in the picture of S. Mathematically, the definition is as 
follows: 

Definition 1.1. A set T is said to be a 2 _n -picture of a bounded set S if: 
(i) S C T, and 
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(ii) T C B(S, 2~ n ) = {x£R 2 : \x - s\ < 2~ n for some s G S}. 

Definition 1.1 means that T is a 2~ n -approximation of S with respect to the Hausdorff metric, 
given by 

d H (S, T) := inf{r : S C B(T, r) and T C B(S, r)}. 

Suppose we are trying to generate a picture of a set S using a union of round pixels of radius 
2~ n with centers at all the points of the form , with i and j integers. In order to draw 

the picture, we have to decide for each pair whether to draw the pixel centered at -^j 
or not. We want to draw the pixel if it intersects S and to omit it if some neighborhood of the 
pixel does not intersect S. Formally, we want to compute a function 

f 1, B((i/2 n J/2 n ),2- n )nS^(& 
f s (n,i/2 n ,j/2 n ) = I 0, fl((i/2»,j/2»),2-2-»)nS = (1.1) 

( or 1, in all other cases 

The time complexity of S is defined as follows. 

Definition 1.2. A bounded set S is said to be computable in time t(n) if there is a function 
f(n, •) satisfying (1.1) which runs in time t(n). We say that S is poly-time computable if there 
is a polynomial p, such that S is computable in time p(n). 

Computability of sets in bounded space is defined in a similar manner. There, the amount 
of memory the machine is allowed to use is restricted. 

To see why this is the "right" definition, suppose we are trying to draw a set S on a 
computer screen which has a 1000 x 1000 pixel resolution. A 2~ n -zoomed in picture of S has 
0(2 2n ) pixels of size 2 _n , and thus would take time 0(t(n) ■ 2 2n ) to compute. This quantity 
is exponential in n, even if t(n) is bounded by a polynomial. But we are drawing S on a 
finite-resolution display, and we will only need to draw 1000 • 1000 = 10 6 pixels. Hence the 
running time would be O(10 6 • t(n)) = 0(t(n)). This running time is polynomial in n if and 
only if t(n) is polynomial. Hence t(n) reflects the 'true' cost of zooming in. 

1.4 Background in complex analysis 

To make the paper self-contained, we list here a few results from complex analysis that will be 
used later. We refer to [Ahl, Dur, Pom] for a more comprehensive discussion. 

Let Q, C C be a simply-connected planar domain with w £ Q. The Riemann Uniformization 
Theorem states that there is unique conformal map of onto the unit disk D with ip(w) = 0, 
vp'(w) > 0. The number r(£l,w) = \/ip'{w) is called the conformal radius of Q. Roughly 
speaking, r(Q,w) measures the size of £1 as viewed from w: 

Proposition 1.2 (Koebe's Theorem). In this notation we have 

r($l,w) 
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We note the following basic monotonicity property of the conformal radius: 

Lemma 1.3. If£l\ C 0,2, w G then r(Q,\,w) < r(^2,w). 

By a theorem of Caratheodory (see e.g. [Pom]), if the boundary <9il is a Jordan curve, then 
the map tft can be extended to a homeomorphism between the closure of f2 and the closed unit 
disk cl(D). 

Let z* = \j~z be the inversion of z with respect to the unit circle \\z\ = 1}. We will make 
use of the following particular case of the Reflection Principle: 

Lemma 1.4. If J C {\z\ = 1} is an open arc, and (ft is a continuous map on D U J which is 
analytic on D, and (ft(J) C {\z\ = 1}, then the map <1? defined by 



is analytic on the domain DU{|z| > 1} U J. 

In particular, if (ft is a conformal map of D onto a domain Q C D with Jordan boundary, and 
K is an open arc, K C d£l n {\z\ = 1}, then $ is a conformal map of B U {\z\ > 1} U J, where 
J = (ft~~ l (K) ((ft is extendable to cl(D) by Caratheodory theorem). 

Let now be a domain with the boundary d£l consisting of finitely many Jordan curves. 
Let / be a continuous function on dfl. A function u : cl(Q) — > C is a solution for the Dirichlet 
problem with the boundary data /, if 

• u is continuous in cl(fi), 

• u is harmonic on Q (An = d xx u + d yy u = 0), and 

• u(z) = f(z) for z 6 dQ.. 

For any / such a solution exists and is unique. Moreover, there exists a unique family of 
measures uj w ,n on dQ such that for any / G C(dQ), 



The measure uj w ^q, is called a harmonic measure. If one fixes K C dCl, the function w i— > 
^uj,n(-^) is harmonic in Q. 

If i7 is simply-connected, then for a set K C dft, we have 



where ift is the Riemann map of Q onto B with ip(w) = 0. 

The Dirichlet problem can be solved probabilistically. Namely, for w G CI, let B w (t) be the 
two-dimensional Brownian motion started at w, and the exit time be defined by 




(1.2) 




Uw,n( K ) = ^length(V(i-Q), 



T = inf{t : B w (t) & fi}. 
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Then the solution of the Dirichlet problem is given by the following formula of Kakutani (see 
e.g. [GM]): 

u(w)=E[f(B w (T))]. 
Note that the harmonic measure for a set K C d£l is now given by 

oj w ,n(K) = P[B w (T)eK}. 

We will make use of the Maximum Principle for harmonic functions (see [Ahl] ) : 

Lemma 1.5. If u\(z) and U2(z) are two functions which are harmonic in Q, continuous on 
the whole cl(O), and u\{z) > U2(z) for z € dfl, then u\(z) > u 2 (z) for all zgfl. 

An easy consequence is the Monotonicity Property of the harmonic measure (see [Pom]): 

Corollary 1.6. // w £ Qi C Vi 2 , K C dVti n dCl 2 , then 

We will also make use of a Distortion Theorem for conformal maps (see [Dur]): 
Theorem 1.7. If <fi is conformal in the disk {z : \z — w\ < r}, then 

I^MI , r l l r W \\2 ^ W*) - *MI ^ I^H L r2| f~ H h2 (i-3) 

(r + |z — (r — |z — w\Y 

and 

r 2 |^V)l /~, |Z ~"!L < \<t>\*)\ < r 2 W(w)\ , r+ } Z ~ W \ (1.4) 
Irv 71 (r + ^ ~ H) 1 v ;I (r - |z - w\f y ' 

1.5 Results 

In Section 2 we propose a new algorithm for computing the Riemann map. We use the random 
walks solution to the general Dirichlet problem to produce a solution to the uniformization 
problem. This gives an extremely space-efficient algorithm. 

The formulation of the theorem will depend on how the boundary of the uniformized domain 
f2 is specified for our algorithm. Since the domain Q we consider is computable, there exists 
a Turing machine M(n) which for a given n computes a function (1.1). Our algorithm may 
then query M(n) for different values of (i/2 n , j/2 n ) to ascertain whether this particular dyadic 
rational point lies within one-pixel distance from dVt. A formal way of saying this is that our 
algorithm will have an access to an oracle for a function given by (1.1). 

Theorem 1.8. There is an algorithm A that computes the uniformizing map in the following 
sense. 

Let Q be a bounded simply- connected domain, and wo 6 Q. dQ is provided to A by an 
oracle representing it in the sense of equation (1.1). Then A computes the absolute values 
of the uniformizing map 4> '■ (^>^o) — ► 0^>O) with precision 2~ n in space bounded by C • n 2 , 
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and time 2°^ n \ where C depends only on the diameter of £1 and d(wo,dQ). Furthermore, the 
algorithm computes the value of(f>(w) with precision 2~ n as long as \4>(w)\ < l — 2~ n . Moreover, 
A queries dfl with precision of at most 2~°( n \ 

In particular, if dVt is polynomial space computable in space n a for some constant a > 1 and 
time T(n) < 2°( na \ then A can be used to compute the uniformizing map in space (7. n max ( a - 2 ) 
and time 2°^. 

In the scale where the entire boundary is given to us explicitly, and not by an oracle for it, 
we have the following. 

Theorem 1.9. There is an algorithm A' that computes the uniformizing map in the following 
sense. 

Let Q be a bounded simply- connected domain, and wq £ Q. Suppose that for some n = 2 k , 
dQ is given to A' with precision - by 0(n 2 ) pixels. Then A' computes the absolute values of the 
uniformizing map (p : (fi, wq) — > (0,0) within an error of 0(1/ n) in randomized space bounded 
by 0(k) and time polynomial in n = 2 k (that is, by a BPL(n) -machine) . Furthermore, the 
algorithm computes the value of (j>(w) with precision 1/n as long as \<f)(w)\ < 1 — 1/n. 

In Section 3, we show that even if the domain we are uniformizing is very simple com- 
putationally, the complexity of the uniformization can be quite high. Moreover, it might be 
difficult to compute the conformal radius of the domain. 

More specifically, the following theorems are established in the Section 3. 

Theorem 1.10. Suppose there is an algorithm A that given a simply- connected domain Q with 
a linear-time computable boundary and an inner radius > \ and a number n computes the first 
20n digits of the conformal radius r(Q,0), then we can use one call to A to solve any instance 
of a #SAT(n) with a linear time overhead. 

In other words, #P is poly-time reducible to computing the conformal radius of a set. 

Theorem 1.11. Consider the problem of computing the conformal radius of a simply- connected 
domain 0, where the boundary ofQ. is given with precision ^ by an explicit collection ofO(n 2 ) 
pixels. 

Denote the problem of computing the conformal radius with precision by CONF(n,n c ). 
Then MAJ n is AC reducible to CONF(n,n c ) for any < c < 1/2. 

1.6 Comparison with known results. 

The first constructive proof of the Riemann Uniformization Theorem is due to Koebe [Koebe], 
and dates to the early 1900's. Formal proofs of the constructive nature of the Theorem which 
follow Koebe's argument under various computability conditions on the boundary of the do- 
main are numerous in the literature (see e.g. [Cheng, BB, Zhou, Hert]). In particular, Zhou 
[Zhou] and Hertling [Hert] give constructive proofs under computability conditions on the 
boundary similar to those used by us. The question of complexity bounds on the construction 
was raised, in particular, in most of the works quoted above. However, the only result known 
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to us was announced by Chou in [Chou]. He states that in the case when the boundary is 
poly(n) computable, the problem of computation of the mapping is in EXPSPACE(n). 

From the practical (that is, applied) point of view, the most computationally efficient 
algorithm used nowadays to calculate the conformal map is the "Zipper" , invented by Marshall 
(see [Mar] ) . The effectiveness of this algorithm was recently studied by Marshall and Rohde in 
[MR]. The "Zipper", however, falls beyond the theoretical upper bound on the complexity of 
this problem, which we establish in Section 2: in the settings of the Theorem 1.8, it computes 
the uniformizing map in space 2°( na ^ and time 2°( na \ and thus belongs to the complexity class 
EXP. It is reasonable to expect then, that an algorithm can be found in the class PSPACE 
which is more practically efficient than "Zipper" . 

2 Computing the uniformization in polynomial space 

Let f2 be a bounded simply-connected planar domain, let K C £1 be a fixed compact set with 
smooth boundary with dist(if, dVt) > 10 • 2~ n . First we discuss a probabilistic algorithm for 
solving the Dirichlet problem in the domain \ K with precision 2~ n . 

2.1 General Dirichlet problem 

The discrete analogue of the Dirichlet problem can be defined as follows. For H C hi? (h>0), 
the interior of H is defined by Int(H) = {a £ H : a ± h, a ± ih £ H}. The boundary of H is 
defined by dH = hi 2 \ (lnt(#) U lnt{KI? \ H)) . We say that a function u defined on H C hi 2 
is discrete harmonic if for any a £ lnt(H) we have 

u(a) = l/4(u(a + h) + u(a — h) + u(a + ih) + u(a — ih)). 

Let be the standard Random Walk (cf [Spi]) on hi 2 started at w £ H, where H is closed 
(dH C H). Let the exit time N be defined as N = min{n : B™ H} — 1. Let / be a function 
on dH. It is almost obvious that the function 

u(w) = E[f(B%)} 

is discrete harmonic on H. This function is called the solution of the Dirichlet problem with 
boundary data f (cf. the continuous case discussion in subsection 1.4). 

Let be a domain with boundary dfl and / £ C(dQ). For h > define H^ = Qnhl 2 . For 
w £ dHh, let fh(w) = f(z), where z is one of the points on d£l closest to w. The solution Uh 
of the corresponding discrete Dirichlet problem, is called the /i-discrete solution of the original 
continuous Dirichlet problem. 

We need the following easy case of the approximating property of the /i-discrete solutions 
(see [Spi], [Laa]). 

Lemma 2.1. Let Q be a domain. Let f be a continuous and locally constant function on dQ, 
taking only and 1 values, and u is the solution of the corresponding Dirichlet problem. Let 
h < 10~ 3 be such that for any z\,Z2 £ dtl with \z\ — z%\ < \fh we have f(z±) = f{z2). Let Uh 
be the h-discrete solution. Then if dist(w, d£l) > \fh, then \u(w) — Uh(w)\ < 2y/h. 
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Since the exit probabilities of a random walk can be computed by a BPL(/i _1 ) machine; 
if the values of / and the boundary dH^ are given by an oracle, then u can be computed in 
the randomized space 0{— log ft,) and time 0(ft~ 2 ). Thus Lemma 2.1 immediately implies the 
following statement about the solution of the general Dirichlet problem: 

Lemma 2.2. There is a randomized algorithm D that computes a solution of the Dirichlet 
problem in the following sense. 

Let Q be a bounded planar domain and K C £1 be a fixed compact set with smooth boundary 
and dist(K, dQ) > 10-2~ n . Suppose that f is the function which is equal to on d£l and 1 on K. 
Then D computes the solution of the corresponding Dirichlet problem with precision 2~ n , 2~ n - 
away from <9f2 U K in space O(n), and time 2°( n \ The computation is done probabilistically, 
and outputs the correct value within an error of 2~ n with probability > \ . 

In particular, if both K and d£l are computable in space n a for some constant a > 1 and 
time T(n) < 2° (na \ Then we can compute the solution of the Dirichlet problem for any point, 
which is at least 2~ n away from dQ. and K in space 0(n a ), and time 2°^T(n). 

2.2 The conformal radius 

Let wq € fl, and let tp be the conformal mapping of Q onto the unit disk B with ^(wq) = 
and i/}'(wq) > 0. Assume that <9S1 is given to us within an error of 2~ n in Hausdorff metric and 
that d(wo, dO,) > 1. As a first application of Lemma 2.2 let us give an algorithm for calculating 
|V>'(wo)| with precision 2~ n in space 0{n a ), and time 2°^T{n). Denote 

wi = w + e~ n and K x = B(w , e~ 2n ) 

Lemma 2.3. Let hi be the solution of the following Dirichlet problem: 



B(0, e~ 2 7(l + e- 2n )V (wo)) C ^{K^ C B(0, e' 2n /(l + e^V (wo)) (2.1) 

Let Bx = V -1 (B{0,e- 2n {l-3e- 2n )ip'{wo))) and B 2 = ip' 1 (5(0, e~ 2n (l - 3e- 2n )^'(w ))) . 
Since 1/(1 + e~ 2n ) 2 < (1 - 3e~ 2n ) and (1 + 3e~ 2n ) < 1/(1 - e" 2n ) 2 , (2.1) implies 



h\{w) = 1, \w — wq\ = 
< hi(w) =0, w £ dQ 
> Afti(w)=0, ro£l]\Jl'i 



Then 




Proof. By the first statement of Theorem 1.7, 



B 1 CK 1 C B 2 



(2.2) 



The functions 



log | ip {w)\ 



-In + log(l - 3e" 2n ) + log ip' (w ) 
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and 



2n + log(l + 3e" 2n ) + logtp'(w ) 

are harmonic in 0, \ B\ and 0, \ B2 respectively, equal to on d£l, and equal to 1 on the 
boundaries of B\ and B2 respectively. By the Maximum Principle, Hi < hi < H2, or, more 
explicitly, 

l0g ^ < hi(w) < - log| y' (2.3) 



-2n + log(l - 3e" 2n ) + log ip'(w ) ~ ~ -2n + log(l + 3e~ 2n ) + log ip'(w ) 

Another application of the same distortion theorem yields 

e" n (l - 3e- n )/(u> ) < \ip(wi)\ < e" n (l + 3e" n )^' \w ) . (2.4) 

Evaluating both sides of the inequality (2.3) at the point wi using (2.4) completes the proof 
of the lemma. □ 

It now follows from Lemma 2.2 that we can compute 4>'{wq) with the same complexity 
constraints as in Lemma 2.2. 

2.3 The Riemann map 

Let hi, Ki be as in the previous section. 

Lemma 2.4. Let \w — wq\ > e~ n . Then 

|log|V>H| - hi(w)(logip'(w ) - 2n)\ < 3 • e~ 2n . 

Proof. By the equation (2.3), 

/ii(w)log(l -3e" 2 ") < log|V»H| - hi(w)(logip'(w ) - 2n) < hi{w) log(l + ?,e~ 2n ) 

To prove the lemma it suffices to notice that hi{w) < 1 and | log(l + x)\ < \x\. □ 

Using Lemmas 2.3 and 2.2, we see that 1^(^)1 is computable with the same restrictions as in 
Lemma 2.2, provided that dist(w, dQ) > e~ n and \w — wq\ > e~ n . 

Now we have to compute axg(i/j(w)) . To achieve this, we introduce another Dirichlet 
problem. Let K2 = B(wq + e~ 2n , e _4ri ), and let h 2 be the solution of the following Dirichlet 
problem: 

h 2 (w) = 1, \w — wq — e~ 2n \ = e 
< h 2 {w) =0, w £ dtt 
K Ah 2 {w)=0, wen\K 2 

Define 

~ {w) _ VH - e~ 2n i;'(wo) 



-An 



ip(w)e 2n ip'(wo) 
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Figure 1: Computing &rg(ip(w)) 



tp is also a Riemann map from Q onto D. Let w 2 '■= ip (0) = ?/> 1 (e 2n ^'(wo))- By Distortion 
Theorem 1.7, 



|u> 2 -u> -e- 2n | < 2-e" 4n . 



As in Lemma 2.4, 



log I i> (w) I - h 2 (w) (log t/r (tu 2 ) - 4n) 



< 3 • e 



-in 



Another application of the Distortion Theorem 1.7 yields | log \ ip'(w2) \ — log \ tp'(wo + e 
2e" 4n . So, finally, 



-2nN 



< 



log \i>(w)\ - h 2 {w)(\og \tp'(w + e 2n )| - 4n) 



< 5 • e 



-'If! 



Now we use a standard formula from hyperbolic trigonometry (see [Thu]) 

cosh C — cosh A cosh B 



where 



cosarg?/>(w) 



A = log- B = l OJ 

1 - h/uw) 



sinh j4 sinh B 

1 + e - 2n ^'(o) 



1 - e - 2n V'(0) 



, and C = log 



1 + \tp(w) 



1 



-2n 



sinhj4 > e n when \w\ > e n , cosh 5 — 1 ~ e 2n . Using the error 



See Figure 1. 

Note that sinh B ' 

estimate in the Lemma 2.4, we obtain that the formula allows us to compute cos arg ip(w) up 
to e~ n , provided that 1^(^)1 < 1 — e _n . 

Using the same argument for computing cos&Tg(ifj(w)/i), we can completely determine the 
value of arg </>(«;). 

Now we can give an algorithm which satisfies the conditions of Theorems 1.8 and 1.9. 
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Figure 2: A a 



Proof of Theorems 1.8 and 1.9. We can create a poly(n) x poly(n) matrix M representing the 
transition probabilities between the poly(n) possible states of the random walk. Simulating 
the random walk for t = poly(n) steps amounts to approximating M*. The required precision 
is also inverse polynomial in n. By Theorem 1.1, this can be done in time polynomial in n, 
and space <9(log 2 n), which imply Theorem 1.9. By changing the scale, and replacing n with 
2™, we obtain Theorem 1.8. □ 



3 Lower bounds on the complexity of uniformization 

In this section we establish Theorems 1.10 and 1.11. 

Let us first remark that by Distortion Theorem 1.7, any algorithm computing values of the 
uniformization map will also compute the conformal radius with the same precision. 

Let A a be the domain D\{|z — 1| < a} - the unit disk with a small bump of radius a 
removed (see Figure 2). 

Fix a large n G N. Let now for < I < 2 n , and let = e 27 ™'/ 2 ™ A 2 -ion be the rotated domain 
A 2 -ion. For a set L = {h, l 2 , . . . , h] with all < h < l 2 ■ ■ ■ h < 2™, let £l L = Q h n Q.i 2 n . . . Qi k . 
Thus VLl is the unit disk with k relatively "spread out" bumps removed. 

Theorem 3.1. For large enough n, 

|r(fi L ,0) - l + A;2- 20n - 1 | < ^ 2 ~ 2 ° n - 

To prove Theorem 3.1, we estimate the conformal radius of A a for an arbitrary a. 

2 -2a 

Lemma 3.2. The conformal radius of A a is equal to 



2 - 2a + a 2 ' 
As a consequence we get that for large n, 

|r(A 2 -ion) - 1 + 2- 20n ~ 1 \ < 2- 30n + 2 (3.1) 

Proof of Lemma 3.2. Let P = C\{Imz = 0, Rez < 0} be the complex plane with the negative 
real axis removed. 
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The function 

maps D conformally onto P , x(0) = 1- It a l so maps A a onto A' b = P \ {\z\ < b}, where 

\ 2 

a 



2 - a, 
Observe also that 

z + b 2 /z - 2b 



h(z) 



(i - by 



maps A' b conformally onto P, with h(l) = 1, h'(l) = ——7 = ^° ~*~ - ■ 

Thus the map <f>o(z) = x 1 h -1 o x(z) maps ID conformally onto A a , and the conformal 
radius of A a is equal to 

□ 

For a set L, let 0l be the conformal map of D onto VLl with </>l(0) = 0, 4>' L (0) > (</>l 
is the inverse of the uniformization map). Let V = (hih,- ■ ■ ,h) be the set L with the first 
element removed. Let g(z) = <p~^} o 4>l(z) be the conformal map of D onto 

r = D\ 4>- L } ({\z\ < 1, \z- e 2rf i/ 2n | < 2- 10n }) . 
Let us also introduce two domains 

r+ = D\{z : \w - z\ < 2~ 10n (l - 2~ 2n )} 

and 

T_=B\{z : \w - z\ < 2- Wn {l + 2- 2n )}, 

where w = (f>j}(l). 

We will use the following property of T 

Lemma 3.3. 

r_crc r+. 

Let us first show how to derive Theorem 3.1 from Lemma 3.3. 

By Lemma 3.3 and Lemma 1.3 (monotonicity of conformal radius), 

r(r_)< 5 '(o)=r(r)<r(r + ) 

Now Lemma 3.2 implies that for large n 

| r (F_) - 1 + 2- 20 "- 1 ! < 2 ~ 22n+2 , |r(r+) - 1 + 2~ 20n - 1 | < 2" 22n+2 , 
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Figure 3: The map (fry and domain T 



and thus 

| 5 '(0) _i + 2- 20n - 1 | < 2~ 22n+2 . 

Note now that 4>l(z) = (fry ° 9{ z )i so r {^L) = g'(0)r(O,y). The Theorem easily follows from 

this relation by induction on the size of L. 

So to establish Theorem 3.1, it is enough to prove Lemma 3.3. 

Proof of Lemma 3.3. Without loss of generality we can assume that l\ = 0. 

Let T = Bfl {\z - 2" 7n+1 | < 1 - 2" 7n }. Note that T C Qy, since V . Let ifr be the 
conformal map of B onto T with tfr(0) = 0, ifr'(0) > 0. 

Let K be the arc [1, e ni2 ~ 7 "}. Note that K C dTndn L >. Let K' = (frj}(K). The normalized 
length of K', \K'\, is the harmonic measure of K in Qy evaluated at zero. By monotonicity 
of harmonic measure (Lemma 1.6) it is bounded above by 2~ 7 ™~ 1 = \K\ and below by the 
harmonic measure of K in T evaluated at zero. 

So 

2 -7n-i > > > 2' 7 "- 1 (l - 2~ 5n ) (3.2) 

The same estimate applies to K" = (fr^i ( [e~ m2 ?n , 1] ) . 
Using these estimates we see that the arc 

j e -7ri2- 7n (l-2- 5 ") e 7ri2- 7n (l-2- 5n )j-j 
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Figure 4: <^y in the neighborhood of w 



is mapped by 4>l' inside K'UK" C dO. Thus we can use the Reflection Principle 1.4 to extend 
the map to a map G of the whole disk {\z — w\ < 2 _7n ~ 1 (l - 2~ 5n )} with G(w) = 1. 

Let J e be the arc [e~ Je , e le ]. Using the fact that T C 0, L i and the monotonicity of harmonic 
measure, we obtain 

\G-\J £ )\ = |^(J e )| = > ^o,r(J e ) = |V-~V0I > 1^1(1 - 2" 5n ) (3.3) 

Letting e — > 0, we obtain 

Now we can use the Distortion Theorem 1.7 applied to the disk {\z-w\< 2- 7n (l - 2~ 5n )} 
to see that 

G({z : \w-z\ < 2- 10 "(l - 2- 2n )}) c{|z-l| < 2~ Wn } 

and 

{|z- 1| < 2- 10n } C : \w-z\ < 2- 10n (l + 2' 2 ")}) 

But this is precisely the statement of the lemma. □ 

Now we are in the position to prove Theorems 1.10 and 1.11. 

Proof of Theorem 1.10. For a propositional formula $ with n variables, let L C {0, 1, . . . , 2 n — 
1} be the set of numbers corresponding to its satisfying instances. Then the boundary of 
Q.L is computable in linear time, given the access to <£. Theorem 3.1 now implies that using 
r(Q,L, 0) we can evaluate \L\ = k, and solve the #SAT problem on <E>, which is exactly Theorem 
1.10. □ 

Proof of Theorem 1.11. Suppose that we are given a string s of n = 2 k zeros and ones. We 
can view it setLc{0,l,...,2 fe -l}. tt L can be obtained from L by a trivial one-layered 
circuit with just NOT gates. Theorem 3.1 implies that using r(O,L,0) with 2~°^ precision, 
we can evaluate \L\ and solve the MAJ„ problem on s, which is exactly Theorem 1.11. □ 
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